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Abstract 

The production of forward jets has been measured in deep inelastic ep collisions at HERA. 
The results are presented in terms of single differential cross sections as a function of the 
Bjorken scaling variable (xsj) and as triple differential cross sections d 3 a /dxBjdQ 2 dp 2 - et , 
where Q 2 is the four momentum transfer squared and p\ - et is the squared transverse mo- 
mentum of the forward jet. Also cross sections for events with a di-jet system in addition 
to the forward jet are measured as a function of the rapidity separation between the forward 
jet and the two additional jets. The measurements are compared with next-to-leading order 
QCD calculations and with the predictions of various QCD-based models. 
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1 Introduction 



The hadronic final state in deep inelastic scattering offers a rich field of research for QCD phe- 
nomena. This includes studies of hard parton emissions which result in well defined jets, pertur- 
bative effects responsible for multiple gluon emissions and the non-perturbative hadronisation 
process. 

HERA has extended the available region in the Bjorken scaling variable, xsj, down to val- 
ues of XBj — 10~ 4 , for values of the four momentum transfer squared, Q 2 , larger than a few 
GeV 2 , where perturbative calculations in QCD are expected to be valid. At these low xsj val- 
ues, a parton in the proton can induce a QCD cascade, consisting of several subsequent parton 
emissions, before eventually an interaction with the virtual photon takes place (Fig. [l]). QCD 
calculations based on "direct" interactions between a point-like photon and a parton from an 
evolution chain, as given by the DGLAP scheme [1-5], are successful in reproducing the strong 
rise of F 2 (xBj, Q 2 ) with decreasing x B j over a large Q 2 range [6-9]. The DGLAP evolution 
resums leading log(Q 2 ) terms. This approximation, however, may become inadequate for small 
XBj, where \og(l/x) terms become important in the evolution equation. In this region the BFKL 
scheme [10-12] is expected to describe the data better, since this evolution equation sums up 
terms in log(l/x). 

Significant deviations from the simple leading order (LO) DGLAP approach are observed 
in the fractional rate of di-jet events [13-15], inclusive jet production [16, 17], transverse energy 
flow [18, 19] and p t spectra of charged particles [20]. Extending the calculations from LO to 
next-to-leading order (NLO) accuracy accounts for some of the deviations observed in jet pro- 
duction, but at low XBj and low Q 2 the description of the measurements is still unsatisfactory. 
Next-to-next-to-leading order (NNLO) calculations do not exist so far and therefore higher or- 
der contributions can only be approximated by phenomenological QCD models, based on LO 
matrix element calculations together with parton shower evolution. Ascribing partonic structure 
to the virtual photon and thus considering so called resolved photon processes, including parton 
showers from both the photon and the proton side, results in an improved description of the data 
including particle production in the forward region (the angular region close to the proton beam 
direction) [13,21-26]. The colour dipole model (CDM) [27,28], which assumes gluon emis- 
sions to originate from independently radiating colour dipoles, is in fairly good agreement with 
the measurements. This suggests that different parton dynamics, not included in the DGLAP 
approximation, are responsible for the observed deviations. 

The large phase space available at low x B j makes the production of forward jets a particu- 
larly interesting process for the study of parton dynamics, since jets emitted close to the proton 
direction lie well away in rapidity from the photon end of the evolution ladder (Fig.[l]). Here a 
new measurement of forward jet production is presented using data collected in 1997 with the 
HI detector, comprising an integrated luminosity of 13.7 pb _1 . The enlarged statistics allows 
to study more differential distributions than previously presented [29-31]. The proton energy is 
820 GeV and the positron energy is 27.6 GeV which correspond to a centre-of-mass-energy of 
v/i « 300 GeV. 

Measurements are presented in regions of phase space where the DGLAP approximation 
might be insufficient to describe the parton dynamics. In inclusive forward jet production this 
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Figure 1 : Schematic diagram of ep scattering with a forward jet taking a fraction Xj et = 
Ej e t/E p of the proton momentum. The evolution in the longitudinal momentum fraction, x, 
from large Xj et to small X gj IS indicated. 



is expected to be the case when the transverse momentum squared of the jet and the photon vir- 
tually are of similar order. More exclusive final states, like those containing a di-jet system in 
addition to the forward jet (called '2+forward jet' ), provide a further handle to control the parton 
dynamics. The forward jet measurements are compared to LO and NLO di- and three-jet cal- 
culations, and different phenomenological QCD models. This measurement is complementary 
to a similar measurement of 7r°-production in the forward direction, which has been presented 
in [26]. 



The conventional description of the parton cascade is given by the DGLAP evolution equations. 
The basic assumption is that the leading contribution comes from cascades with strong ordering 
in the virtualities of the parton propagators in the evolution chain, with the largest virtualities 
reached in the hard scattering with the photon. This implies strong ordering of the transverse 
momenta of the emitted partons (k t ). Since their virtualities and transverse momenta squared are 
small compared to the hard scale, Q 2 , the propagators can be treated as massless and assumed 
to be collinear with the incoming proton (collinear approach). The interaction is assumed to 
take place with a point- like photon (DGLAP direct). 

If the scale of the hard subprocess is larger than Q 2 , the structure of the virtual photon might 
be resolved and the interaction take place with one of the partons in the photon. In this case a 



2 QCD Models and Theoretical Calculations 
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partonic structure is assigned to the photon and a photon parton density function is convoluted 
with the matrix element, which within the DGLAP model means that two evolution ladders are 
introduced, one from the photon side and one from the proton side of the hard subprocess. This 
is called the resolved photon model (DGLAP resolved) and is described in [22, 23]. 

The BFKL ansatz predicts strong ordering in the longitudinal momentum fraction of the 
parton propagators but no ordering in their virtualities. This means that the virtualities and 
the transverse momenta of the propagators can take any kinematically allowed value at each 
splitting. One consequence of this is that the matrix element must be taken off mass-shell and 
convoluted with parton distributions which take the transverse momenta of the propagators into 
account (unintegrated parton densities). 

The CCFM equation [32-35] provides a bridge between the DGLAP and BFKL descriptions 
by resumming both \og(Q 2 ) and log(l/a;) terms in the relevant limits, and is expected to be 
valid in a wider x range. The CCFM equation leads to parton emissions ordered in angle. An 
unintegrated gluon density is used as input to calculations based on this model. 

A different approach to the parton evolution is given by the colour dipole model (CDM), 
in which the emissions are generated by colour dipoles, spanned between the partons in the 
cascade. Since the dipoles radiate independently there is no ordering in the transverse momenta 
of the emissions and the behaviour of the parton showers is in that sense similar to that in the 
BFKL case. 

The measurements performed here are compared to several QCD models: 

• The RAPGAP [36] Monte Carlo program, which uses LO matrix elements supplemented 
with initial and final state parton showers generated according to the DGLAP evolution 
scheme for the description of DIS processes (RG-DIR). It can be interfaced to HERA- 
CLES [37], which simulates QED-radiative effects. RAPGAP also offers the possibility 
to include contributions from processes with resolved transverse virtual photons (RG- 
DIR+RES). 

• The DJANGO [38] program with the CDM as implemented in ARIADNE [39]. Param- 
eters of ARIADNE are tuned using the CTEQ6M [40] parton density functions and the 
data sets [19,20,41]. 

• The CASCADE Monte Carlo program [42, 43], which is based on the CCFM formal- 
ism [32-35]. Two different versions of the unintegrated gluon density are used, J 2003- 
set-1 and set-2 [44]. The difference between these two sets is that in set-1 only singular 
terms are included in the splitting function, whereas set-2 also takes the non-singular 
terms into account. These unintegrated gluon densities are determined from fits to the 
F 2 (x, Q 2 ) data obtained by HI and ZEUS in 1994 and 1996/97. 

Simulated events from the RAPGAP (RG-DIR) and DJANGO Monte Carlo programs are 
processed through the detailed HI detector simulation [45] in order to test the understanding of 
the detector and to extract correction factors. 

The forward jet cross sections are compared to LO (a s ) and NLO (a 2 ) calculations of di-jet 
production via direct photon interactions as obtained from the DISENT program [46,47] . Since 



6 



the jet search is performed in the Breit frame the selected events always contain at least one jet 
in addition to the forward jet, such that comparisons with the DISENT predictions are ade- 
quate. The renormalisation scale (fj%) is given by the average p\ of the di-jets from the matrix 
element ((^di-jets))' wni l e the factorisation scale {jA) is given by the average pi of all forward 
jets in the selected sample 1 j et )). The calculations are corrected for hadronisation effects, 
which are estimated using CASCADE together with the KMR parton density function [48]. The 
KMR parton density function takes only the matrix element and one additional emission into 
account and should therefore be suitable for correcting the NLO di-jet calculations. The correc- 
tion factors for hadronisation effects (1 + <5 H ad) are determined by calculating the ratio bin- wise 
between the hadron and parton level cross sections, obtained using the same jet algorithm and 
kinematic restrictions. 

In the analysis of events with two jets in addition to the forward jet, the measured cross 
sections are compared to the predictions of NLOJET++ [49]. This program provides pertur- 
bative calculations of cross sections for three-jet production in DIS at NLO (cuf) accuracy. In 
NLOJET++, where the factorisation scale can be defined for each event, fj% and ^ are set to 
the average p\ of the forward jet and the two hardest jets in the event. The NLOJET++ cal- 
culations are corrected to hadron level using CASCADE together with the unintegrated gluon 
density J2003 set-2 [44]. 

The NLO calculations by DISENT and NLOJET++ are performed using the CTEQ6M [40] 
parametrisation of the parton distributions in the proton. The uncertainty in the NLO calcu- 
lations originating from the PDF uncertainty is estimated by using the CTEQ eigenvector sets 
according to [40]. The scale uncertainty for these calculations is estimated by simultaneously 
changing the renormalisation and factorisation scales (fj%, pjj) by a factor of 4 up and 1/4 down. 
In CASCADE the renormalisation scale is changed by the same factors and in each case 
the unintegrated gluon density is adjusted such that the prediction of CASCADE describes the 
inclusive F 2 data [50, 51]. The forward jet cross section is then calculated to estimate the upper 
and lower limit of the scale uncertainty. The resulting uncertainty in the cross section prediction 
is less than 10% at the smallest xsj and decreases for higher XBj (these errors are not shown 
in the figures). The parton densities and the scales used in the QCD calculations are given in 
tabled] 





CASCADE 


RG-DIPv/RES 


DISENT 


NLOJET++ 


Mr 




Q' 2 + (Ptdi-iets) 


(Ptdi-iets) 


(Ptjetl + Pt,jet2 + Pt,fwdjet)/3 




s + Q 2 


Q 2 + (pmhJ 


(Phet) 


kPtjetl + Pt,jet2 + Pi.fwdjei)/^ 


proton PDF 


J2003 set-1 k-2 


CTEQ6L [40] 


CTEQ6M 


CTEQ6M 


photon PDF 




SaSID [52] (RES onfy) 







Table 1: The renormalisation and factorisation (pyj) scales, and the parton density func- 
tions used in the different programs. The average squared transverse momentum of the forward 
jet, (pfj et ), is 45 GeV 2 for the single differential forward jet cross section, and 24, 55 and 183 
GeV 2 for the three different pf -bins in the triple differential cross sections, s is the invariant 
mass squared of the di-quark system. 

'For the triple differential forward jet cross section, d 3 a/dxBjdQ 2 dp 2 ,- et , this means different factorisation 
scales for the three different ptjet bins. 
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In [53] next-to-leading order calculations of the forward jet cross section are presented, in 
which the contributions from direct and resolved virtual photons are taken into account in a 
consistent way. The inclusion of NLO contributions from the resolved part corresponds to an 
additional gluon emission in a direct process and thus may constitute an approximation of the 
NNLO direct cross section. 



3 The HI Detector 

A detailed description of the HI detector can be found in [54-56]. The detector elements 
important for this analysis are described below. The coordinate system of HI is defined such 
that the positive z axis is in the direction of the incident proton beam. The polar angles are 
defined with respect to the proton beam direction. 

The interaction vertex is determined with the central tracking detector consisting of two 
concentric drift chambers (CJC) and two concentric z drift chambers (CIZ and COZ). The 
kinematic variables x and Q 2 are determined from a measurement of the scattered electron in the 
lead- scintillating fibre calorimeter (SpaCal) and the backward drift chamber (BDC), covering 
the polar angular range 153° < 6 < 177°. 

The SpaCal has an electromagnetic section with an energy resolution of 7%/ ' sj ' E /GeV © 
1%, which together with a hadronic section represents a total of two interaction lengths. Iden- 
tification of the scattered electron is improved using the BDC, situated in front of the SpaCal. 
The scattering angle of the electron is determined from the measured impact position in the 
BDC and the reconstructed primary interaction vertex. 

The hadronic final state is reconstructed with the Liquid Argon calorimeter (LAr), the cen- 
tral tracking detector and the SpaCal. The LAr calorimeter is of a sandwich type with liquid 
argon as the active material. It covers the range 4° < 9 < 154°. In test beam measurements 
pion induced hadronic energies were reconstructed with a resolution of about 50% / a/ E/GeV © 
2% [57]. The measurement of charged particle momenta provided by the central tracking de- 
tector is performed in a solenoidal magnetic field of 1.15 T with a precision of a p /p 2 = 0.003 
GeV -1 . 

The luminosity is determined from the rate of Bethe-Heitler events (e + p^e + 7 + p) with 
a precision of 1.5%. 

The scattered electron is triggered by its energy deposition in the SpaCal. For events used 
in this analysis, with the electron energy required to be above 10 GeV, the trigger efficiency is 
essentially 100%. 



4 Experimental Strategy and Phase Space Definition 

Differences between the various approaches to the modelling of the parton cascade dynamics are 
most prominent in the region close to the proton remnant direction, i.e. away from the photon 
side of the ladder. This can be understood from the fact that the strong ordering in virtuality of 
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the DGLAP description gives the softest /^-emissions closest to the proton whereas in the BFKL 
model the emissions can be arbitrarily hard in this region, as long as they are kinematically 
allowed. 

In most of the HERA kinematic range the DGLAP approximation is valid. A method to 
suppress contributions from DGLAP like events is to select events with a jet close to the proton 
direction (a forward jet) with the additional constraint that the squared transverse momentum 
of this jet, Pt t j et , is approximately equal to the virtuality of the photon propagator, Q 2 (see 
Fig. [l])- This will suppress contributions with strong ordering in virtuality as is the case in 
DGLAP evolution. If, at the same time, the forward jet is required to take a large fraction of 
the proton momentum, Xj et = Ej et /E p , such that Xj et ^> XBj, the phase space for an evolution 
with ordering in the longitudinal momentum fraction, as described by BFKL, is favoured. By 
requiring Xj et ^> XBj contributions from zeroth order processes are also suppressed. Based on 
calculations in the leading log approximation of the BFKL kernel, the cross section for DIS 
events at low XBj and large Q 2 with a forward jet [58, 59] is expected to rise more rapidly with 
decreasing XBj than expected from DGLAP based calculations. 

DIS events are selected by requiring a scattered electron in the backward SpaCal calorimeter 
and a matching track in the backward drift chamber (BDC), applying the following cuts: 

E' e > 10 GeV 
156° <6 e < 175° 
0.1 < y < 0.7 
0.0001 < x Bj < 0.004 
5 GeV 2 < Q 2 < 85 GeV 2 

where E' e is the energy of the scattered electron, 6 e its polar angle, and y is the inelasticity vari- 
able. The lower cut on Q 2 and and the upper on y reduce the background from photoproduction. 

Jets are defined using the fc^-jet algorithm [60,61] with combined calorimeter and track 
information [62] as input (applied in the Breit-frame). The selection further requires the recon- 
struction of at least one jet in the laboratory frame, satisfying the cuts below: 

p tJet > 3.5 GeV 
7° < 9 jet < 20° 
x jet > 0.035 



where the ptjet- and jet -cuts are applied in the laboratory frame. If there is more than one jet 
fulfilling these requirements the most forward is chosen. For the single differential cross section 
measurement an additional cut 0.5 < r = p 2 j et /Q 2 < 5 is applied. 

Data are presented as single differential cross-sections as a function of x B j (da/dxBj), and 
triple differential cross-sections as a function of xbj in bins of Q 2 and p 2 j et 
(d 3 a/dxBjdQ 2 dp 2 j et ). Another event sample, called the '2+forward jet' sample, is selected 
by requiring that, in addition to the forward jet, at least two more jets are found. Out of these, 
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the two with the highest transverse momenta are chosen. This provides further constraints on 
the kinematics at the expense of reducing the data sample. 

For the '2+forward jet' sample the p t is required to be larger than 6 GeV for all 3 jets. 
The other cuts on the forward jet are kept the same as specified above, and no p 2 j et / Q 2 -c\xi 
is applied. The two additional jets are required to lie in pseudorapidity, r] = — lntan(6 l /2), 
between the electron and the forward jet, r] e < i]j et i < r]j et2 < r]f wd j et . 

The final numbers of events used for the single and the triple differential forward jet cross 
section are 17316 and 23992, respectively. The number of selected '2+forward jet' events is 
854. 



5 Correction Factors and Systematic Uncertainties 

The RAPGAP and DJANGO programs, together with a simulation of the HI detector, are used 
to correct the data for acceptances, inefficiencies, and bin to bin migrations due to the finite 
detector resolutions. The shapes of the distributions of the DIS kinematic variables and the jet 
variables for the forward jet sample, as defined in section HJ are compared to the predictions 
from RAPGAP and DJANGO. This is done by reweighting the Monte Carlo XBj distributions 
to give the best possible agreement with data and by studying how well the distributions of 
the other kinematic variables are described. The distributions are reproduced equally well by 
the predictions of RAPGAP and DJANGO after the detector simulation. In Fig. detector level 
distributions are shown for x Bj , Ej et and p\j et / Q 2 for the forward jet samples, with and without 
the p 2 

jet/ Q 2 -cut applied. These distributions are normalised to the number of events and thus 
give a shape comparison to investigate the understanding of the detector, independently of the 
normalisation of the physics models. 

The hadron level cross sections are extracted by applying correction factors to the data in 
order to take detector effects into account. The correction factors are calculated as the ratio of 
the CDM Monte Carlo prediction at the hadron and detector levels, in a bin-by-bin procedure. 
These factors correct the data from detector level to non-radiative hadron level, i.e. the data are 
also corrected for QED radiative effects. RAPGAP and CDM give similar values over the full 
kinematic range covered in this investigation. The correction factors are generally between 0.7 
and 1.2 but in a few kinematic bins they reach values of 0.5 or 1.4 due to limited resolution of 
the jet quantities. The variations in the correction factors between the two Monte Carlo models 
are included in the systematic error. 

The purity and acceptance 2 are found to be larger than 30% in all bins. For the '2+forward 
jet' analysis they are larger than 40% in all bins. 

The systematic errors are estimated for each data point separately as the quadratic sum of 
the individual errors described below. The following systematic errors are considered: 

2 The purity (acceptance) is obtained from the same Monte Carlo simulations as used for the correction factors 
and is defined as the number of simulated events which originate from a bin and are reconstructed in it divided by 
the number of reconstructed (generated) events in that bin. 
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Forward jets 





Figure 2: Control plots for the forward jet selection. The sample with no p 2 j et /Q 2 -cut applied 
(upper) and the sample with the 0.5 < p 2 j et /Q 2 < 5-cut applied (lower) are shown. The 
distributions are at detector level and normalised to unity. All variables are measured in the 
laboratory frame. Comparisons are made to the predictions of the DJANGO (full line) and 
RAPGAP ( dashed line ) Monte Carlo programs. 



• The hadronic energy scale uncertainty is determined to be 4%. In order to estimate the 
related uncertainty of the measured forward jet cross section, the reconstructed hadronic 
energies in the DJANGO/ ARIADNE simulation were increased and decreased by this 
amount. The average resulting error is typically 8% for both the single differential for- 
ward jet cross section and the triple differential forward jet cross section, and 13% for the 
'2+forward jet' cross section. 

• The electromagnetic energy scale as measured in the SpaCal is known to an accuracy of 
1%. Changing the scale by this amount in the forward jet cross section calculations using 
DJANGO/ ARIADNE results in an average systematic error of typically 3% for the single 
and triple differential measurement, and 1% for the '2+forward jet' measurement. 

• The uncertainty on the measured scattering angle of the electron is estimated to be 
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1 mrad, which contributes typically 1% to the error in the forward and '2+forward jet' 
cross section. 

• The error from the model dependence is taken as the difference between the correction 
factors calculated from the DJANGO/ ARIADNE and the RG-DIR Monte Carlo pro- 
grams. Taking this variation into account yields a systematic error of about 5% for the 
single differential forward jet cross section, 8% for the triple differential case and 13% in 
the '2+forward jet' cross section. 

• The PHOJET [63, 64] Monte Carlo generator was used in order to estimate the extent 
to which DIS forward jet events could be faked by photoproduction (Q 2 ~ GeV 2 ) 
background. The influence on the measurement is found to be negligible. The error 
attributed to this source of uncertainty is taken to be 1%. 

• The uncertainty of the luminosity measurement is estimated to be 1.5%. 

The averages of these sums are 10%, 12% and 14% for the single differential, triple differ- 
ential and the '2+forward jet' cross section, respectively. In the figures the systematic errors 
due to the energy scale uncertainty of the calorimeters (A Systl ) are shown separately as bands 
around the data points, whereas the other systematic errors (A Syst2 ) are included in the error bars 
together with the statistical errors. The errors are given separately in the tables. 

6 Results 

6.1 Single Differential Cross Section 

The measurement of the single differential forward jet cross section is presented at the hadron 
level in the phase space region defined in section |4j The phase space for DGLAP evolution is 
suppressed by the additional requirement 0.5 < p 2 j et /Q 2 < 5 as discussed in section|4| 

The measured single differential forward jet cross sections are listed in table |2j In Fig. [3^. 
they are compared with LO (a s ) and NLO (a 2 ) calculations from DISENT. The calculations 
are multiplied by (1 + <5had) to correct to the hadron level. The uncertainty from the factori- 
sation and renormalisation scales, and the uncertainty in the PDF parametrisation, are added in 
quadrature to give the total theoretical error, which is shown as a band around the histogram 
presenting the theoretical prediction. In Fig. and c the data are compared to the various QCD 
models. 

In Fig.|3k it can be observed that, at small XBj, the NLO di-jet calculations from DISENT 
are significantly larger than the LO contribution. This reflects the fact that the contribution 
from forward jets in the LO scenario is suppressed by kinematics. For small XBj the NLO 
contribution is an order of magnitude larger than the LO contribution. The NLO contribution 
opens up the phase space for forward jets and improves the description of the data considerably. 
However, the NLO di-jet predictions are still a factor of 2 below the data at low XBj- The 
somewhat improved agreement at higher xbj can be understood from the fact that the range in 
the longitudinal momentum fraction which is available for higher order emissions decreases. 
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Figure 3: The hadron level cross section for forward jet production as a function ofxsj com- 
pared to NLO predictions from DISENT (a) and to QCD Monte Carlo models (b and c). The 
shaded band around the data points shows the error from the uncertainties in the calorimetric 
energy scales. The hatched band around the NLO calculations illustrates the theoretical uncer- 
tainties in the calculations, estimated as described in the text. The dashed line in (a) shows the 
LO contribution. 



From Fig.|3b it is seen that the CCFM model (both set-1 and set-2) predicts a somewhat 
harder x Bj distribution, which results in a comparatively poor description of the data. 

Fig- Et shows that the DGLAP model with direct photon interactions alone (RG-DIR) gives 
results similar to the NLO di-jet calculations and falls below the data, particularly in the low 
XBj region. The description of the data by the DGLAP model is significantly improved if contri- 
butions from resolved virtual photon interactions are included (RG-DIR+RES). However, there 
is still a discrepancy in the lowest £ Bj -bin, where a possible BFKL signal would be expected 
to show up most prominently. The CDM model, which gives emissions that are non-ordered in 
transverse momentum, shows a behaviour similar to the RG DIR+RES model. Analytic calcu- 
lations where resolved photon contributions are included to NLO order [53] again give similar 
agreement with the data as the RG DIR+RES model [22]. 



6.2 Triple Differential Cross Sections 

In this section data are presented as triple differential forward jet cross sections. The total 
forward jet event sample is subdivided into bins of Q 2 and p 2 j et - The triple differential cross 
section da j dxBjdQ 2 dp 2 j et versus xsj is shown in Figs. |4||6] for three regions in Q 2 and p 2 j et . 
Fig. ID presents the cross section compared to NLO (a 2 ) calculations, including theoretical er- 
rors, represented by error bands. In Fig. 13 and |6] comparisons to QCD Monte Carlo models 
are shown. The same parton density functions and scales are used as in the measurement of the 
single differential cross section. The cross section values are listed in table |3l 

From Fig.|4]it can be observed that the NLO calculations in general undershoot the data but 
similarly to the single differential cross section the NLO calculations get closer to the data at 
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higher x B j and so too, due to the kinematics, at higher Q 2 . The NLO calculations also give a 
better description of data for harder forward jets. In the highest • t -bin the difference between 
data and NLO is less than the (large) uncertainty in the NLO calculations in several x B j-bim. 
This is consistent with the results from a previous measurement on inclusive jet production [17]. 
A possible explanation is that jets with high p t remove a large fraction of the energy from the 
parton ladder, leaving limited energy available for additional emissions. Thus, the parton ladder 
is shorter and more like the NLO configuration. The general trend is that NLO calculations 
agree better with data as x B j , Q 2 and p 2 j et increase. For high Ptjet the phase space for LO starts 
to open up, which also makes the NLO prediction more reliable. 

The comparisons between data and QCD models are discussed in three different kinematic 
regions as specified below. These regions are however not strictly separated, but overlap. In 
all three regions the CDM and DGLAP resolved (RG-DIR+RES) models give very similar 
predictions (see Fig.0) indicating that a breaking of the ordering of the virtuality is necessary 
to describe the data. As already observed in the single differential measurement the CCFM 
model predicts a somewhat harder x B j distribution than seen in the data. This is true for the full 
kinematic range and leads to the poor description of the data as seen in Fig. |5] 

Pt, jet ~Q 2 (r~ 1) 

In this region events with parton emissions ordered in p t are suppressed, and thus parton dy- 
namics beyond DGLAP may show up. The data are best described by the DGLAP resolved 
model (RG-DIR+RES) as observed in Fig.HJ) and f. 

Phet <Q 2 (r< 1) 

The region where Q 2 might become larger than p 2 j et is dominated by direct photon interactions. 
However, since r can take values up to 1.8 in the most DGLAP-like bin (Fig.lU), events with 
Ptj et of the same order or even greater than Q 2 are also contributing. This gives an admixture 
of events with emissions non-ordered in virtuality. This may explain why the DGLAP direct 
model (RG-DIR), although closer to the data in this region than in others, does not give good 
agreement with the data except for the highest x B j-b'm. The CDM and DGLAP resolved model 
(RG-DIR+RES) reproduce the data very well in this region. 

Pt, jet >Q 2 (r> 1) 

The kinematic region where p 2 j et is larger than Q 2 is typical for processes where the virtual 
photon is resolved. As expected the DGLAP resolved model (RG-DIR+RES) provides a good 
overall description of the data, again similar to the CDM model. However, it can be noted that 
in the regions where r is the highest and x B j small, CDM shows a tendency to overshoot the 
data. DGLAP direct (RG-DIR) gives cross sections which are too low (see Fig.|6]d, g and h). 

6.3 Events with Reconstructed Di-jets in Addition to the Forward Jet 

Complementary to the analyses reported in sections |6~T1 and |6"21 where the ratio p 2 j et /Q 2 has 
been used to isolate regions where a possible BFKL signal is enhanced, another method is used 
to control the evolution kinematics in the analysis reported here. By requiring the reconstruction 
of the two hardest jets in the event in addition to the forward jet, different kinematic regions can 
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be investigated by applying cuts on the jet momenta and their rapidity separation as described 
in more detail in section HI 

In this scenario it is demanded that all jets have transverse momenta larger than 6 GeV. By 
applying the same p t j et -cut to all three jets, evolution with strong A; r ordering is not favoured. 
Decreasing the p t j et -cut is not possible in this analysis due to detector resolutions. The jets are 
ordered in rapidity according to i]f wd j et > r)j et2 > r\jet x > r] e with rj e being the rapidity of the 
scattered electron. The cross section is measured in two intervals of Ar/i = rjj et2 — r)j etl . If 
the di-jet system originates from the quarks q± and q 2 (see Fig.|7]), the phase space for evolution 
in x between the di-jet system and the forward jet is increased by requiring that A771 is small 
and that Ai] 2 = i]f wd j et — r]j et2 is large. A771 < 1 favours small invariant masses of the di-jet 
system and thereby small values of x g (see Fig. |7]). With Ar] 2 large, x g carries only a small 
fraction of the total propagating momentum, leaving the rest for additional radiation. It should 
be kept in mind, however, that only the forward jet is explicitely restricted in rapidity space, by 
the demand that it has to be close to the proton axis. The directions of the other jets are related 
to the forward jet through the Ar] requirements. When Ar] 2 is small, it is therefore possible that 
one or both of the additional jets originate from gluon radiation close in rapidity space to the 
forward jet. With Ar/i large, BFKL-like evolution may then occur between the two jets from 
the di-jet system, or, with both Ar/i and Ai] 2 small, even between the di-jet system and the 
hard scattering vertex. By studying the cross section for different At] values one can test theory 
and models for event topologies where the k t ordering is broken at varying locations along the 
evolution chain. 

The cross sections for events containing a di-jet system in addition to the forward jet are pre- 
sented as a function of Ar/ 2 in Figs.lSlfTUIfor all '2+forward jet' events , and for the requirements 
Ar]i < 1 and Ar/i > 1, respectively. The measured cross sections are given in table HJ For the 
Ar]i < 1 region the cross section falls at low Ar/ 2 since the phase space becomes smaller when 
the 3 jets are forced to be close together. Fig.[8]gives a comparison of data to NLO (a^) predic- 
tions with theoretical error contributions included as bands. In Figs. |9] and [10] comparisons to 
QCD models are presented. 

In this investigation the same settings of the QCD models are used as in sections 16.11 and 
16.21 while the NLO three-jet cross sections are calculated using NLOJET++. 

From Fig. [S] it is observed that NLOJET++ gives good agreement with the data if the two 
additional hard jets are emitted in the central region (Ar] 2 large). It is interesting to note that 
a fixed order calculation (af), including the log (l/x) -term to the first order in a s , is able to 
describe these data well. However, the more the additional hard jets are shifted to the forward 
region (Ai] 2 small), the less well are the data described by NLOJET++. This can be understood 
from the fact that the more forward the additional jets go, the higher the probability is that one 
of them, or even both, do not actually originate from quarks but from additional radiated gluons. 
For gluon induced processes, which dominate at small x, NLOJET++ calculates the NLO con- 
tribution to final states containing one gluon jet and two jets from the di-quarks, i.e. it accounts 
for the emission of one gluon in addition to the three jets. Thus, events where two of the three 
selected jets originate from gluons are produced by NLOJET++ only in the real emission cor- 
rections to the three-jet final state, which effectively means that these kinematic configurations 
are only produced to leading order (af). The most extreme case, where all three reconstructed 
jets are produced by gluons, is not considered by NLOJET++. This results in a depletion of the 
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theoretical cross section in the small Af] 2 region, which is more pronounced when Ar}% is also 
small, i.e. when all three jets are in the forward region. Consequently a significant deviation 
between data and NLOJET++ can be observed for such events (see the lowest bin in Fig.|Sj5). 
Accounting for still higher orders in a s might improve the description of the data in this domain 
since virtual corrections to the production of two gluons could increase the cross section for 
such final states, and additional gluon emissions would enhance the probability that one of the 
soft radiated gluons produces a jet that fulfills the transverse momentum requirement applied in 
this analysis. 

For the '2+forward jet' sample CCFM is not describing well the shape of the redistributions 
(Fig.|9k, b and c). 

As explained above, evolution with strong A; t -ordering is disfavoured in this study. Radiation 
that is non-ordered in k t may occur at different locations along the evolution chain, depending 
on the values of Ar/i and Ar/ 2 . As can be seen from Fig. 10, the colour dipole model gives good 
agreement in all cases, whereas the DGLAP models give cross sections that are too low except 
when both A771 and Ar/ 2 are large. For this last topology all models and the NLO calculation 
agree with the data, indicating that the available phase space is exhausted and that little freedom 
is left for dynamical variations. 

If one or both jets from the di-jet system are produced by gluon radiation, which is in- 
creasingly probable the more forward these jets go, it necessarily means that the k t ordering 
is broken. In this context it is noteworthy that CDM provides the best description of the data 
while the other models, including the DGLAP-resolved model, fail in most of the bins. The 
'2+forward jet' sample differentiates CDM and the DGLAP-resolved model, in contrast to the 
more inclusive samples where CDM and RG-DIR+RES give the same predictions. The conclu- 
sion is that additional breaking of the k t ordering is needed compared to what is included in the 
resolved photon model. 



7 Summary 

An investigation of DIS events containing a jet in the forward direction is presented. Various 
constraints are applied, which suppress contributions to the parton evolution described by the 
DGLAP equations and enhance the sensitivity to other parton dynamics. Several observables 
involving forward jet events are studied and compared to the predictions of NLO calculations 
and QCD models. 

Leading order (a s ) calculations of the single differential forward jet cross section, da/dxBj, 
are well below the measurements, which is expected since forward jet production is kinemati- 
cally suppressed in LO. NLO di-jet calculations improve the description of the data but remain 
too low at small values of xbj- This is also the case for predictions based on the DGLAP di- 
rect model. The DGLAP resolved photon model (RG-DIR+RES) and the colour dipole model 
(CDM) come closest to the data. 

The total forward jet sample is subdivided into bins of Q 2 and pL et such that kinematic 
regions are defined in which the effects of different evolution dynamics are enhanced. In the 
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most DGLAP enhanced region, (Q 2 3> P$j et )> and in the region where contributions from 
resolved processes are expected to become important (p 2 - t ^> Q 2 ), the measured triple differ- 
ential forward jet cross sections are well described by the CDM and the DGLAP resolved model 
(RG-DIR+RES). In the BFKL region (Q 2 ~ p 2 tjet ) the CDM and DGLAP resolved model (RG- 
DIR+RES) again reproduce the data best. A general observeration is that the DGLAP resolved 
model and CDM tend to fall below the data at low x Bj , Q 2 and p 2 . The cross sections predicted 
by the DGLAP direct model (RG-DIR) are consistently too low in all regions and especially at 
low x Bj . 

The NLO di-jet calculations from DISENT describe the data for the largest values of x Bj at 
high values of Q 2 and pf, but fail for low values of these variables. 

The measured cross section for events with a reconstructed di-jet system in addition to the 
forward jet are in good agreement with the predictions of NLOJET++ if the additional jets are 
emitted in the central region. As expected, deviations are observed for other jet topologies. The 
data are best described by the CDM. The DGLAP resolved model (RG-DIR+RES) is below 
the data as is, to an even greater extent, the DGLAP direct model (RG-DIR). This result gives 
the first evidence for parton dynamics in which there is additional breaking of the /c r ordering 
compared to that provided by the resolved photon model. 

The CCFM model, as implemented in CASCADE, with two different parametrisations of 
the unintegrated gluon density, fails to describe the shape of both the single and triple differ- 
ential cross sections, as well as the '2+forward jet' cross section. This might be caused by the 
parametrisation of the unintegrated gluon density and/or the missing contributions from split- 
tings into quark pairs. 

The observations made here demonstrate that an accurate description of the radiation pat- 
tern at small x Bj requires the introduction of terms beyond those included in the DGLAP direct 
approximation (RG-DIR). Higher order parton emissions with breaking of the transverse mo- 
mentum ordering contribute noticeably to the cross section. Calculations which include such 
processes, such as CDM and the resolved photon model, provide a better description of the 
data. The similar behaviour of CDM and the DGLAP resolved model (RG-DIR+RES), which 
describe the data best, indicates that the inclusive forward jet measurements do not give a signif- 
icant separation of the models. However, in the more exclusive measurement of '2+forward jet' 
events a clear differentiation of the models is obtained since, in contrast to CDM, the DGLAP 
resolved model (RG-DIR+RES) fails to describe the data. 
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X B j 


do 1 dx Bj (nb) 


Astat (nb) 


Asysti (nb) 


As ys t2 (nb) 


0.0001-0.0005 


925 


± 17 


+110 

-100 


+77 
-77 


0.0005-0.001 


541 


± 12 


+54 
-55 


+23 
-24 


0.001-0.0015 


264 


1 o 

± 8 


+30 
-28 


+11 
-11 


0.0015-0.002 


153 


± 6 


+19 
-16 


+8 
-8 


0.002-0.003 


74.5 


±3.0 


+10.7 
-8.0 


+1.9 
-1.8 


0.003-0.004 


36.7 


±2.0 


+2.1 
-5.7 


+2.4 
-2.4 



Table 2: Single differential forward jet cross sections in bins ofx Bj- The statistical error (A Stat ), 
the error from the uncertainty of the calorimetric energy scales (As yst i) and from the other 
systematic errors (A S y St2 ) are specified. 
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Figure 4: The hadron level triple differential cross section for forward jet production as a func- 
tion ofxsj, in bins of Q 2 (GeV 2 ) and p 2 j et (GeV 2 ). The data are compared to the prediction 
of NLO (full line) and LO (dashed line) calculations from DISENT. Both calculations are cor- 
rected for hadronisation effects. The band around the data points illustrates the error due to the 
uncertainties in the calorimetric energy scales. The band around the NLO calculations illus- 
trates the theoretical uncertainties in the calculations. In each bin the range in and the average 
value ofr = p 2 j et /Q 2 is shown. 
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Figure 5: The hadron level triple differential cross section for forward jet production as a 
function ofx Bj, in bins ofQ 2 (GeV 2 ) and p\j et (GeV 2 ). The data are compared to the predictions 
of CASCADE. The band around the data points illustrates the error due to the uncertainties in 
the calorimetric energy scales. In each bin the range in and the average value of r = p 2 j et /Q 2 
is shown. 
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Figure 6: The hadron level triple differential cross section for forward jet production as a 
function ofx Bj, in bins ofQ 2 (GeV 2 ) and p 2 j et ( GeV 2 ). The data are compared to the prediction 
of RAP GAP DIR, RAPGAP DIR+RES and CDM. The band around the data points illustrates 
the error due to the uncertainties in the calorimetric energy scales. In each bin the range in and 
the average value ofr = Ptj et /Q 2 is shown. 
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Figure 7: A schematic diagram of an event giving a forward jet and two additional hard jets. 
These may stem from the quarks (qi and q 2 ) in the hard scattering vertex or from gluons in the 
parton ladder. x g is the longitudinal momentum fraction carried by the gluon, connecting to the 
hard di-jet system (in this case qi and g 2 J • 
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Figure 8: The cross section for events with a reconstructed high transverse momentum di-jet 
system and a forward jet as a function of the rapidity separation between the forward jet and 
the most forward-going additional jet, Ar] 2 . Results are shown for the full sample and for two 
ranges of the separation between the two additional jets, Arji < 1 and Ai]i > 1. The data are 
compared to the predictions of a three-jet NLO calculations from NLO JET+ + (1 + 5had)- The 
band around the data points illustrates the error due to the uncertainties in the calorimetric 
energy scales. The band around the NLO calculations illustrates the theoretical uncertainties 
in the calculations. 
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Table 3: Triple differential cross sections in bins of Q 2 , p 2 and XBj- The statistical error 
(A Stat ), the error from the uncertainty of the calorimetric energy scales (A Systl ) and from the 
other systematic errors (A Syst2 ) are specified. 
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Table 4: '2+forward jet' cross sections in bins of Ai] 2 for all Arji, Arji < 1 and Arji > 1. The 
statistical error (A Stat ), the error from the uncertainty of the calorimetric energy scales (A Systl ) 
and from the other systematic errors (A Syst2 ) are specified. 
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